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SUMMARY 

An approximate analysis was developed to analyze the postbuckling behavior 
of through-width delaminations in a laminated coupon. The analysis contains two 
parameters which are determined using a finite element analysis. After calcula- 
ting the parameters for a few configurations, the approximate analysis was used 
to analyze many other configurations. Lateral deflections and mode I strain- 
energy release rates obtained with the approximate analysis were compared with 
results from the finite element analysis. For the configurations analyzed, the 
approximate analysis agreed very well with the finite element results. 

INTRODUCTION 

In laminated composite structures under compression loads, delaminations 
often precipitate failure. Even small, seemingly benign delaminations may induce 
localized buckling, which causes high interlaminar stresses. A buckled delami- 
nated region may grow rapidly and lead to structural instability. To assess the 
criticality of a delamination, an analysis is needed to predict the rate of 
instability-related delamination growth. A key component of any such analysis 
is an accurate geometrical ly-nonl inear stress analysis. Unfortunately, general 
purpose analyses, such as geometrical ly-nonl inear finite element analysis, tend 
to be expensive. Specialized inexpensive approximate analyses are needed. 

This paper presents an approximate analysis for postbuckling of through- 
width delaminations in a laminated composite coupon (fig. 1). This configuration 


-7 XI 





is perhaps the simplest that exhibits instability-related delamination 
growth. Hence, it is an ideal candidate for initial study. Furthermore, 
this configuration was studied in reference 1 using a geometrical ly- 
nonlinear finite element analysis. Thus, reference solutions are available. 

In the following sections the development of an approximate analysis is 
outlined. Lateral deflections and mode I strain-energy release rates obtained 
with the approximate analysis are compared with results obtained with the 
analysis described in reference 1. 

NOMENCLATURE 

2 

cross-sectional area of region A, m 
2a length of delamination before loading, m 

2a' effective length of delamination before loading, m 

2a axial length of delamination after loading, m 

b width of coupon, m 

C^,C 2 constants used in calculating Gj 

E Young's modulus, GPa 

E^,Ejj extensional moduli for regions A and D, GPa 

2 

Gj mode I strain-energy release rate, J/m 

Gp transverse shear modulus of region D, GPa 

2 

G^^ strain-energy release rate due to a moment, J/m 

4 

I moment of inertia, m 

Ip moment of inertia of region D, m^ 

M moment, N*m 

crack closing moment, N*m 
crack opening moment, N*m 


2 



n 




x,y 

6 

/V 

6 

6 




shape factor equal to 1.2 (ref. 2) 

axial loads in regions A, B, C, and D, N 

buckling load calculated using finite element analysis, N 
applied load, N 

applied load corresponding to incipient buckling, N 

thickness of delaminated region, m 

rectangular Cartesian coordinates, m 

maximum lateral deflection of buckled region, m 

lateral deflection corresponding to peak Gj, m 

lateral deflection corresponding to closing of crack tip, m 

axial strain in region A 

midplane axial strain in region D 


ANALYSIS 

Development of Governing Equations 

In this section approximate governing equations are derived for a laminate 
with a postbuckled through-width delamination. The laminate was subdivided 
into four regions, as shown in figure 2. Because of symmetry, only half of 
the laminate was modeled. The laminate was assumed to be of width b. 

Regions B and C are assumed to be perfectly bonded. Regions A and D 
are unbonded. Regions A, B, and C have constant axial strain. Hence, 
the force-displacement relations are those for a simple rod subjected to 
axial load. Region D is assumed to have zero slope at both ends. To 
describe the nonlinear behavior of region D, equations (1) and (2) for post- 
buckling of a column were used. 
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a - a = a 
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ref. (3) 


( 2 ) 


where 6, a, "a, and Pp are lateral deflection, axial length before and after 
deformation, and load, respectively. Equations (1) and (2) were derived using 
a strength of materials analysis of a column. 

To combine regions A, B, C, and D, equilibrium and compatibility 
conditions must be considered. The equilibrium condition for the axial force 
is 


Pa " Pd = 




(3) 


Compatability requires the shortening of regions A and D to be identical. 
Hence, 


Pa ^ 

^A^A 


= a - a 


(4) 


Equations (2), (3), and (4) can now be combined to obtain the governing 
equation for the laminate in terms of one unknown, 6. 
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Equation (5) can be solved explicitly to obtain 6. 



( 6 ) 


In a later section, lateral deflections calculated with equation (6) will be 
compared with results obtained with a finite element analysis. 

Calculation of Mode I Strain-Energy Release Rate 
In reference 1 the mode I strain-energy release rate (Gj) was shown to 
dominate instability- related delamination growth. The complexity of the load 
transfer at the crack tip prevents a simple, exact calculation of Gj. An 
approximate procedure is presented here. 

Transverse normal stress at the crack tip is the net result of two 
opposing processes. The lateral deflection, 6, causes a moment Mg which 
tends to open the crack tip. M^ can be calculated by considering moment 
equilibrium of region D (see sketch below). 


^0 



5 



I 


Mq is found to be 



( 7 ) 


After the delaminated region buckles, an increase in applied load causes no 
change in the load carried by the buckled column (region D in the sketch), 
(see eq. (1)). 


But the load carried by region C continues to increase with increased applied 
load. Hence, load must be transferred from C to A. The eccentricity in the 
load path causes a moment, M^, which tends to close the crack tip. Because 
the closing moment arises from the difference in axial forces in regions C 
and D, is expected to be strongly dependent on the quantity (P^ - Pp). 

Hence, the closing moment is assumed to be 

= -'=1 I Cc - Pd> (8) 

The constant is difficult to determine analytically because it 
depends on the complex stress diffusion process at the delamination front. 
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In the present method, it is left as a free constant to be determined through 
finite element analysis. 

The strain-energy release rate for a cantilever beam loaded by a moment 

O 

is M /2EIb (ref. 4). Assuming region D to behave like a cantilever beam, 
the strain-energy release rate associated with the moments and is 



(M. + M )' 

-w 


( 9 ) 


However, Gj is not necessarily equal to Gj^^. Because the specimen is not 
symmetric about the delamination, the moments can create both normal and 
shear stresses. Assuming the behavior near the crack tip is linear, Gj is 
a constant fraction of G^^. 


G 


I 


c m 


(10) 


Combining equations (6) through (10) and using the rule of mixtures to 
calculate Pq, Gj can be expressed in terms of 6 as 
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LEpAp + E^A^ V^D^D 


(EqAp + EflAfl) + 


A"A^ 


A c 2^,2 

16 a^ 





( 11 ) 


The constants Cj and C 2 in equations (8), (10) and (11) have to be 
calculated with a finite element analysis; the procedure for calculating them 
will be discussed in a later section. 
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Extremum Behavior of Gj 

Reference 1 showed that Gj initially increases with increasing load 
and lateral deflection and then decreases. Since finite elements were used, 
no closed form expression was obtained to describe this behavior. In this 
section equations are presented which describe when Gj first becomes 
nonzero (i.e., initial buckling), the load and deflection for peak Gj, the peak 
value of Gj, and the load and deflection for closure of the crack tip (which 
results in zero Gj). 

Initial buckling occurs when the load carried by the delaminated region, 

Pp, reaches the buckling load. Until buckling occurs, the load carried by 
the delaminated region is linearly related to the remote load. 


p = p 


( 12 ) 


Equating equations (1) and (12) yields an expression for the applied load for 
incipient buckling (and hence nonzero Gj). 









3^AdGd / 


(13) 


The maximum value of Gj occurs when the total peeling moment, (M^ + M^), 
reaches a maximum. The maximum value for (M^ + M^) is found by solving 
equation (14) 





Cit 

2 3Pj 


= 0 


(14) 
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The term is calculated from equation (6); the term 



is 


calculated using equation (3) and the rule of mixtures. Solving equation (14), 
one obtains the lateral deflection corresponding to peak Gj, 6. 


7T A^E^ApEpC^t 


The corresponding applied load is obtained by solving equation (6) for Pj, 
that is 





+ A^E^) + 


16a^ 


(16) 


The peak value of Gj is calculated using 6 in equation (11). 

To determine the load and deflection at which the crack tip closes, 
equation (17) is solved. 


Mq + = 0 (17) 

Substituting expressions for and into equation (17) (using equation (3), 
rule of mixtures, and equation (16) to express Pq as a function of 6) 
results in 


3 16a^Pp(ApEp A^E^) 
TT^AAEAApEpCit 


(18) 
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Note that 6 = 26. The corresponding applied load is calculated by using 
6 in equation (16). 


Calculation of Constants, C-j and Cg 

Finite element solutions are required to determine the constants C-| 
and C 2 * The finite element analysis is discussed briefly in a later 
section. 

To calculate C-j, equation (15) was used. Finite element analysis was 
used to determine the lateral deflection, 6, at which Gj is maximized 
for a particular configuration. 

After calculating C^ , C 2 was calculated using equation (10). was 
calculated with equation (9) and Gj was calculated with the finite element 
analysis. Gj^ and Gj were calculated at the lateral deflection corresponding 
to the maximum Gj for a particular configuration. 

C-| and C 2 could have been calculated using any two points, but the 
calculations would have been more tedious. 

Both constants were found to be independent of the delamination length 
and load; C-j was also independent of the thickness of the buckled region. 
Consequently, Gj can be calculated with the approximate analysis for many 
different configurations after studying a few configurations with the finite 
element analysis. 

Modifications for Short, Thick Delaminated Regions 

If the buckled region is short and thick, the analysis described in the 
preceding sections should be modified to account for rotation at the ends of 
the delaminated region. Clamp conditions were assumed to obtain equation (6). 
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Rotation at the ends of the delamination cause the delamination to behave 
like it is longer than it is. That is, the buckling load is lower and the 
lateral deflection is greater than for the case of clamped ends. The 
effective lengthening can be estimated from the relative magnitudes of Pp 
calculated with equation (1) and the buckling load calculated by a finite 
element analysis. If we ignore the small shear correction term 


nTT^E^X 


in equation (1), Pp is inversely proportional to 



Hence, 


the effective length is approximately 



(19) 


Finite Element Analysis 

A two-dimensional, geometrically nonlinear finite element analysis was 
used to obtain rigorous reference solutions. These solutions were used to 
calculate C-j, C2, and 2a' for use in the approximate solution. The 
solutions were also used to determine the accuracy of the approximate 
analysis. 

The finite element analysis is described in detail in reference 1. A 
typical finite element mesh is shown in figure 3. The mesh contained 813 nodes 
and 740 four-node isoparametric elements. Similar meshes were used for other 
delamination lengths (2a) and depths (t). Strain-energy release rates were 
calculated using the crack closure technique reported in reference 5. 
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Description of Specimen Configuration 
The specimen configuration used in the parametric analysis consisted of 
unidirectional graphite/epoxy (regions C and D in fig. 2) bonded to an 
aluminum bar (regions A and B in fig. 2). The graphite/epoxy had a 
thickness t (see fig. 1). The aluminum had a thickness of 6.0 mm. The 
Young's moduli and shear modulus were assumed to be: 

aluminum: E = 67 GPa (9.7 x 10^ PSI) 

graphite/epoxy: E = 140 GPa (20 x lO® PSI) 

G = 5.9 GPA (.85 x 10® PSI) 

RESULTS AND DISCUSSION 

The objective of this section is to illustrate the potential of the 
approximate analysis for analyzing specimens containing through-width 
delaminations. Specimens with different delamination lengths, delamination 
depths (t), and applied loads were analyzed. A finite element analysis was 
used to obtain reference solutions. 

Two calculated parameters were examined — lateral deflections and mode I 
strain-energy release rate (Gj). The lateral deflection was considered 
an important parameter for study because intuitively one might expect the 
severity of the interlaminar stresses to be related to the degree of post- 
buckling. Also the accuracy of the calculated lateral deflection (determined 
by comparison with finite element results) is a measure of how well the 
gross deformation behavior of the specimen is modeled. Gj was considered 
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because results presented in reference 1 suggested that delamination growth 
rates are dominated by Gj. 

In the results presented, corrections for end rotation are not included 
except where specifically indicated. This was done so that the simplest 
form of the approximate analysis could be evaluated. 

Figures 4 through 6 show calculated lateral deflections 6 for various 
delamination lengths, delamination depths, and applied loads. The approximate 
analysis (eq. (6)) and the finite element analysis agree very well for most 
of the cases. The differences are primarily a result of rotation at the 
ends of the delamination. The rotation is greatest for the case t = 0.762 mm. 
Since ignoring end rotation "stiffens" the system, it is not surprising 
that the finite element analysis always predicted a lower buckling load and 
a greater deflection. 

Figure 7 shows that for t = 0.762 mm the approximate and finite element 
analyses agree very well if corrections are included for end rotation. For 
both 2a = 25.4 mm and 2a = 38.1 mm the difference between the effective 
length and actual length of the delamination was calculated (with eq. (19)) 
to be approximately 

(2a‘ - 2a) = 2t (20) 

Figures 8 to 11 compare mode I strain-energy release rates calculated 
with the approximate analysis (eq. (11)) and the finite element analysis. 

Figure 8 shows the relationship between Gj, lateral deflection, and 
delamination length for t = 0.508 mm. The constants C-j and C 2 were 
calculated (as described earlier) using only the finite element results for 
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the case 2a = 25.4 mm. The figure shows the approximate analysis agrees 
very well with the finite element results for all four delamination lengths. 

Figures 9 through 11 show the relationship between Gj, applied load, 
and delamination length. Three values of t are examined: 0.508, 0.254, 
and 0.762 mm. The value of the constant C-j determined from figure 8 (in 
which t = 0.508 mm) was used for all three values of t. C 2 had to be 
calculated for each value of t. The finite element results for the shortest 
delamination length (for each t) were used to calculate C 2 . 

The analyses agree very well for the cases t = 0.508 mm (fig. 9) 
and 0.254 mm (fig. 10). For t = 0.762 mm (fig. 11) the agreement is fair. 
Recall that when lateral deflections were calculated, the approximate 
analysis was not as accurate for the case t = 0.762 mm. However, with the 
corrections for end rotations, the approximate analysis performed very well 
even for t = 0.762 mm. If these corrections are used, the Gj calculations 
also are greatly improved, as shown in figure 12. To determine how these 
corrections might affect calculations for other values of "t," the curves 
in figure 9 were recalculated using the corrections for end rotations. 

Equation (20) was used to determine the effective length (instead of eq. (19)). 
Comparison of figures 9 and 13 shows the approximate analysis agrees more 
closely with the finite element analysis if end rotation effects are included. 

CONCLUDING REMARKS 

An approximate analysis for postbuckling of a through-width delamination 
in a laminated composite coupon was developed. Lateral deflections and 
mode I strain-energy release rates (Gj) obtained with the approximate 


14 



analysis were compared with results from a geometrical ly-nonl inear finite 
element analysis. 

Specimens with different delamination lengths, delamination depths, and 
applied loads were considered. In most cases, lateral deflection and Gj 
obtained with the approximate analysis agreed very well with the finite 
element results. For the cases in which the agreement was only fair, 
excellent agreement was obtained by incorporating corrections for rotation 
of the ends of the delamination. 

Before calculating Gj v/ith the approximate analysis, two constants 
had to be determined from a finite element analysis. Hence, the approximate 
analysis cannot stand alone (except for calculating lateral deflections). 
However, a significant advantage of the approximate analysis is that the 
effect of various parameters on Gj can easily be determined from the 
governing equations. 
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Figure 7.- Lateral deflection vs. load for two delamination lengths. Two versions 
of the "Strength of materials" solution are shown: (a) without correction for end 
rotations (b) with correction for end rotations (t = .762 mm). 
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Figure 8.- Relationship between Gj, lateral defh 

length, (t = . 508 mm) . 
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Figure 13. - Gi versus applied load for sevei 
The "Strength of materials" solution has 
rotation ( t = .508 mm ). 
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